Joint multipartite photon statistics by on/off detection 
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We demonstrate a method to reconstruct the joint photon statistics of two or more modes 
of radiation using on/off photodetection performed at different quantum efficiencies. The 
two-mode case is discussed in details and experimental results are presented for the bipartite 
states obtained after a beam-splitter fed by a single photon state or a thermal state. 

o 

The reconstruction of the joint photon distribution of two or more correlated modes of radiation 
plays a crucial role in fundamental quantum optics 0] and finds relevant applications in quantum 
communication Q], imaging 0] and spectroscopy 0|. Nevertheless, photodetectors suited for this 
' purpose are currently not available, since the few existing examples |5( still suffer from limitations. 

CN ■ On the other hand, reconstruction by quantum tomography [fj] is not an easily implementable 

technique suited for a widespread use. 

Recently 0], a maximum-likelihood (ML) method based on on/off detection performed at dif- 
ferent quantum efficiencies has been developed and demonstrated for reconstructing the photon 
distribution of single-mode states. The results are reliable and accurate also for relatively low 



o 
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\Q \ quantum efficiency of the detector. Since for many applications multipartite states are needed, in 

this letter we extend our previous results to this case as well. In particular, the bipartite case (easily 
extendable to multipartite) will be discussed in details. Examples of experimental reconstructions 
for bipartite state are presented to test and assess our method. 

The statistics of on/off detection performed with quantum efficiency r] on a single-mode state 
p is given by po n = Yln^v^Pn \pin = 1 — Pov] where p n = (n\g\n) is the photon distribution 
(diagonal matrix elements) of the state, and A nn = (1 — r/) n . By performing independent on/off 
photodetection on two (spatially separated) modes of radiation, globally described by the two-mode 
O" 1 density matrix g, the joint on/off statistics is given by 



POOr] = y^A^ArjkQnk (1) 
nk 

POln = y^r;n(l ~ A?fc)gnfc (2) 
nk 

PlOn = y,(l — A vn )A v kQnk , (3) 
nk 

and, of course, p\\ n = 1 — poo?? — PlOn — Poiv where g n k = {{nk\g\nk}} (\nk)) = \n)®\k}) is the joint 
photon distribution of the two modes. Once the value of the quantum efficiency is known, the above 
equations provide a relation between the statistics of clicks and the actual statistics of photons. 
At a first sight this represents a scarce piece of information about the state under investigation. 
However, if the on/off statistics is collected for a suitably large set of efficiency values, then the 
information is enough to reconstruct the joint photon distribution of the bipartite state. We adopt 
the following strategy: by placing in front of the detector K filters with different transmissions, 
we may perform the detection with K different values rj u , v = 1...K, ranging from t\\ — ^/min 
to a maximum value r]x = f?max equal to the nominal quantum efficiency of the detector. Upon 
writing g = (p 00m , p 00m , p 01m , . . . , p 01m ,p 10m , . . . , p WrjK ) and q = ( g 00 , g 01 , g 10 , . . . ) , according to 
Qnk -> Q P with p = 1 + k + n(l + N), i.e. k = (p - 1) mod (1 + N) and n = (p - 1 - k)/(l + N) 
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we can summarize the on/off statistics as 



1,..SK p =!,.,(! +Nf 



(4) 



where we have introduced the matrix B 



Afin-^fik I 1 — 1) •■> K 
B w = { A Mn (l - A„ k ) fj, = K + 1, .., 2K 
(l-A„ n )A^ fi = 2K + l,..,3K 



(5) 



If the Q n kS are negligible for n,k > N and the 77^'s are known, then Eq. represents a finite 
statistical linear model for the positive unknown q p . The maximum-likelihood (ML) solution of 
this LINPOS problem is well approximated by the iterative algorithm [HI] 
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In Eq. (jbj) qy denotes the p-th element of reconstructed statistics at the i-th step, 
9iA{lp }] theoretical on/off probabilities as calculated from Eq. Q at the i-th step, 
whereas hn are the measured frequency of the events with quantum efficiency ry^, i.e h = 

/n^, ra M being the total num- 
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(foOm , ■ ■ ■ , f00ri K , hlrn , ■ ■ ■ , hlr) K , /lOr?i , • • • , fl®r) K ) , with fij Vfl 

ber of runs performed with rj = ry^. 

The convergence of the algorithm may be checked by the total error 



which measures the distance of the reconstructed statistics (of clicks) from the measured one: the 
algorithm is stopped when e reaches its minimum, or goes below a certain threshold value. 




FIG. 1: Schematic diagrams of the experimental setups. On the left: setup to reconstruct the joint photon 
statistics of the bipartite (entangled) states obtained by splitting a PDC-II heralded single-photon state 
by balanced and unbalanced beam-splitters. On the right: setup to reconstruct the joint photon statistics 
of a bipartite (classically correlated) state obtained by splitting a multithcrmal state by a balanced beam 
splitter. 



In order to test our method we have experimentally performed the reconstruction for different 
bipartite states. The schematic diagrams of the experimental setups are reported in Fig. El As a 
first test, we use our method to reconstruct the bipartite (entangled) states obtained at the output 
of a beam splitter fed by the heralded single-photon state generated by parametric down conversion 
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FIG. 2: Reconstruction of the joint photon distribution of a bipartite state: PDC heralded photon state 
split by a beam splitter (balanced and unbalanced). Dashed lines represent one standard deviation. 

(PDC). In our set-up a 0.2 W, 398 nm pulsed (with 200 fs pulses) laser beam, generated by second 
harmonic of a titanium - sapphire beam at 796 nm, pumps a type II BBO 5x5x1 mm crystal. Upon 
detecting a photon in a branch of the degenerate PDC emission one triggers the presence of the 
correlated photon in the other direction. The heralded single-photon is then impinged onto a beam 
splitter (BS) with unexcited second port, thus generating a bipartite entangled state of the form 
\ip)) = i/r|01}) — t] 10)), t being the transmissivity. After the BS both arms are measured by 
an on/off detectors. All the detectors were APD silicon photodetectors, whose quantum efficiency 
has been calibrated with the traditional PDC scheme |l2j. The proper set of quantum efficiencies 
is obtained by inserting before the BS several Schott neutral filters (NF) of different transmittance, 
evaluated by measuring the ratio between the counting rates with the filter inserted and without 
it. The data for the reconstruction have been taken using K — 34 values of tj from T/min — 

0.015 

to 

J/max — 0.325. In correspondence of the detection of a photon in arm 1, a coincidence window 
has been opened on both detectors on arm 2. This is obtained by sending the output of the first 
detector as start to two Time to Amplitude Converters (TAC) that receive the detector signal 
as stop. The 20 ns window is set to not include spurious coincidences with PDC photons of the 
following pulse (the repetition rate of the laser is 70 MHz) . The TAC outputs are then addressed 
both to counters and to an AND logical gate for measuring coincidences between them. These 
outputs, together with one TAC Valid Start (giving us the total number of opened coincidence 
windows), allows to evaluate the frequencies on/off needed for reconstructing the joint photon 
statistics of the bipartite state. The background has been evaluated and subtracted by measuring 
the TAC and AND outputs out of the triggered window. 

In order to verify the method in different cases we considered 4 different alternatives given 
by the combination of a balanced (r = 0.5) or unbalanced (r = 0.4) BS with either large band, 
red glass filters (RG) with cut-off wave length at 750 nm, or interference filters (IF), with peak 
wave length at 796 nm and a 10 nm FWHM. The reconstructed statistics for these four situations 
are shown in Fig. El The uncertainties have been evaluated as described in Ref. '7\. As it is 
apparent from Fig. |2] the reconstructed state well corresponds to a single photon in one of the 
modes. Only the elements £01 ) £10 are different (within uncertainties) from zero and their ratio is 
the value expected by the ratio of the outputs ports of the BS (unity for the balanced one, 2/3 
for the unbalanced one). As expected in this regime, no multi-photon component is observed: i.e. 
£n> 0O2> and so on are zero within uncertainties. The small uncertainties show that also less 
unbalanced BS would be distinguishable. 
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As a second example we consider a single branch of PDC emission without triggering, which 
corresponds to a multithermal state with number of modes of the order of ~ 10 3 . A bipartite state 
is generated by impinging this signal onto a beam splitter with the second port unexcited. The 
output bipartite state is classically correlated (not entangled, but not factorisable) , with the two 
partial traces corresponding to multithermal states. The expected on/off statistics is given by 

Pmn = M M G" + (7) 
Polv = ^[^ + ir rN)-^-(fx + r,N)-^ (8) 
Pior, = ^ [(H + ri(l - t)N)-» -(h + rjN)-L 1 } (9) 

respectively, where N is the average number of photons and \x the number of modes. 




FIG. 3: Reconstruction of the joint photon distribution of a bipartite state: attenuated PDC multithermal 
distribution split by a balanced beam splitter. The two histograms correspond to reconstructed distributions 
(dark gray) compared to multithermal ones. On the bottom left: fidelity F of the reconstructed marginal 
distributions to multithermal ones as a function of the number of the iterations of the algorithm. The inset 
show the total error e. The last box show the measured frequencies /oo, /oij /io as function of the quantum 
efficiency compared with the expected ones (multithermal, solid lines). 

In this case the state has been produced by pumping a 5x5x5 mm type I BBO crystal by a 
beam of a Q-switched triplicated (to 355 nm) Neodimium-Yag laser with pulses of 5 ns, power up 
to 200 mJ per pulse and 10 Hz repetition rate. Due to the very high power of the pump beam, a 
state with a large number of photons is generated. We have therefore attenuated (by 1 nm FWHM 
IF and neutral filters) the multithermal state before splitting and detection. Measurements at the 
different quantum efficiencies has been obtained by inserting (before the BS) Schott neutral filters, 
whose calibration has been obtained by measuring the power of a diode laser (at the same wave 
length of the used PDC emission) before and after them. The data for the reconstruction has been 
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taken using K — 35 values of tj from T/min — 

0.05 to 

??max — 0.25. The coincidence scheme has 
been realized by addressing two Q-switch triggered pulses to two TAC modules as starts, and the 
detectors outputs as stops. Then, having set properly a 20 ns coincidence window, we sent the 
two TAC outputs to an AND logic port, and the Valid Stops to counting modules (together with 
one TAC Valid Start and the AND output). The results of this reconstruction are shown in Fig. 
01 Also in this case the comparison among theoretical expectations and reconstructed statistics 
is rather good. The fidelity F = Y^ n V QnQ™ th of the reconstructed distribution to the expected 
multithermal {g™ th }is larger than 99% for both the marginals. Notice that the optimal number 
of iterations (i.e. leading to maximum average fidelity of the two marginals) corresponds to the 
minimum of e, thus confirming the good convergence properties of the algorithm. 

In conclusion, we demonstrated a method to reconstruct the joint photon statistics of two or 
more modes using on/off photodetection. Experimental reconstruction have been presented for 
the bipartite states obtained after a beam-splitter fed by a single photon state or a thermal state. 
Our results clearly show that the ML reconstruction based on on/off detection can be successfully 
applied to measure the joint photon statistics for multipartite systems. 
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PRIN 2005023443-002 and 2005024254-002), by Regione Piemonte (E14), and by "San Paolo foun- 
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